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Abstract 



Hydrogen and/or helium accreted by a neutron star from a binary companion may undergo thermonuclear fusion. At different mass ac- 
cretion rates different burning regimes are discerned. Theoretical models predict helium fusion to proceed as a thermonuclear runaway 
for accretion rates below the Eddington limit and as stable burning above this limit. Observations, however, place the boundary close 
to 10% of the Eddington limit. We study the effect of rotationally induced transport processes on the stability of helium burning. For 
the first time detailed calculations of thin helium shell burning on neutron stars are performed using a hydrodynamic stellar evolution 
code including rotation and rotationally induced magnetic fields. We find that in most cases the instabilities from the magnetic field 
provide the dominant contribution to the chemical mixing, while Eddington- Sweet circulations become important at high rotation 
rates. As helium is diffused to greater depths, the stability of the burning is increased, such that the critical accretion rate for stable 
helium burning is found to be lower. Combined with a higher heat flux from the crust, as suggested by recent studies, turbulent mixing 
could explain the observed critical accretion rate. Furthermore, close to this boundary we find oscillatory burning, which previous 
studies have linked to mHz QPOs. In models where we continuously lower the heat flux from the crust, the period of the oscillations 
increases by up to several tens of percents, similar to the observed frequency drift, suggesting that this drift could be caused by the 
cooling of deeper layers. 
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1. Introduction 

Neutron stars in low-mass X-ray binaries (LMXBs) accrete mat- 
ter from a companion star through Roche-lobe overflow. Due 
to conservation of angular momentum, the matter forms a disk 
around the neutron star, where it radiates away a large frac- 
tion of the potential energy before reaching the neutron star sur- 
face. The companion is typically a main sequence star, donating 

, matter with a large hydrogen abundance. For a number of sys- 
tems which likely have a white dwarf companion, hydroge n may 
be deficient but helium not (Nelson etal. 1986; Webbin k et al.l 

H983l:lin'tZand et al.ll2005h . 

The accreted matter is thought to quickly spread over the 
neutron star surface, forming a thin layer. When this layer is a 
few meters thick, the temperature and density at the base reach 
the ignition conditions for hydrogen and helium fusion. The ther- 
monuclear burning can proceed in a stable or an unstable man- 
ner, depending on the temperature dependence of the burning 
and cooling rates, which depend on the conditions in the thin 
burning s hell. These conditions are larg ely set by the mass accre- 
tion rate ('Fuiim oto et"aDl 19811: see also iBildstenll 1998 ). Models 
predict that hydrogen burning via the CNO cycle is stable for 
M > 10~^MEdd (assuming spherical symmetry), while it is un- 
stable for lower accretion rates (assuming a neutron star mass of 
1.4M0 with a 10km radius, a hydrogen mass fraction X = 0.7 
and a CNO mass fraction Zcno =0.01; MEdd is the critical ac- 
cretion rate where the accretion luminosity equals the Eddington 
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luminosity). For helium burning via the triple-alpha process, 
theory plac es the critical accretion rate at the E ddington limit 
M ^ MEdd (Fuii moto et aP 1 1981 ': "Ayasli & Joss T982), which 
suggests that stable helium burning does not take place in prac- 
tice. If the accreted matter is hydrogen-deficient, the transition 
is expected to take place at an even higher accretion rate due 
to the lack of h eating by stable hydrogen burning: M lOMEdd 
(lBildstenlll998i) . 

Nevertheless, there are observational indications that stable 
heliu m burn ing takes place at lower accretion rates. Van Paradijs 
et al. (Il988h observe from ten sources, most of which are thought 
to accrete hydrogen-rich material, that the fraction of hydrogen 
and helium that is burned in a stable manner increases strongly 
with the mass accretion rate, w hile no bursts are observed at ac- 
cretion rates of M > O.SMEdd- ICornelisse et al.l (|2003) find for 
4U 1820-30, which likely has a hydrogen-deficient mass donor, 
that the burst rate is lower by a factor of 5 for M > O.OTMEdd 
than for lower mass accretion rates, while no bursts are observed 
at M > 0.14MEdd- Therefore, the transition of unstable to sta- 
ble burning likely takes place within several tens of percents of 
O.lMEdd. 

A possible solution for this discrepancy is the inclusion of 
mixing processes. If hydrogen and/or helium are mixed to larger 
depths, where temperature and density are higher, the conditions 
under which the fuel is burned are different. An efficient source 
of mixing is formed by ro t ationally induced hydrodynamic in- 
stabilities (Fujimoto""l993|). lYoon et al.1 (l2004a) studied the ef- 
fect of rotational mixing processes on thin shell burning in ac- 
creting white dwarfs a nd found burning to proceed more stably. 
iPiro & BildstenI (l2007l) calculated the effect of turbulent mixing 
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in neutron stars , including m ixing due to rotationally induced 
magnetic fields ('Spruit"2002'). They find that the magnetic fields 
are the main driver of the mixing. Their analysis of helium accre- 
tion onto an iron layer shows that in the case of unstable burning, 
rotationally induced mixing allows for lower recurrence times 
for helium flashes. Furthermore, their models including mix- 
ing allow for steady-state burning at accretion rates below MEdd- 
Analytically the dependence of A^st, the accretion rate at the tran- 
sition between stable and unstable burning, on the rotation rate 
of the neutron star, Q., is derived as M^t Qp-^^. For neutron 
stars that rotate slower than 10% of the Keplerian velocity at the 
surface, Piro & Bildsten" ('2007') find 0.25 MEdd < M^t < MEdd- 

Piro & Bildsten ( 2007) created one-dimensional numerical 
models of the neutron star envelope using analytical approxima- 
tions of chemical mixing and of helium burning. In this paper 
we study the effect of mixing due to rotation and a rotationally 
induced magnetic field on the stability of helium burning in the 
neutron star envelope in more detail. We present the results of 
detailed numerical simulations, where we use a one-dimensional 
multi-zone model of the outer layers of the neutron star. We in- 
troduce both rotation and magnetic fields in the model. Using a 
stellar-evolution code we investigate the stability of the burning 
as a function of the mass accretion rate, the rotation rate of the 
neutron star core and the heat flux emanating from the crust. 

2. Method 

To model the outer layers of the neutron star we employ a stel- 
lar evolu tion code that is a modified version of the code used by 
lYoon et al.l (|2004a) to model nuclear burning of helium shells in 
white dwarfs. It implicitly solves the stellar structure equations 
on a one-dimensional grid consisting of a series of mass shells. 
The effect of the centrifugal force on the stellar structure is taken 
into account. Consequently, the mass shells are isobars rather 
than spherical shells (cf. iHeger et al.ll2000l) . The chemical evo- 
lution is modeled by an extensive nuclear network. The effects 
of rotation on the chemical mixing and the transport of angular 
momentum are calculated. The effects of a rotationally induced 
magnetic field on the mixing of an gular mome ntum and chem- 
ical species are calculated as well (iSpruitI 120021) . Convection is 
taken into account using the Ledoux criterio n Neu trino losses 
are calculated using the results of Itoh et al l (Il996h . We make 
use of the OPAL tables (Iglesias & Rogers 1996) for determin- 
ing opacities. Heating due to viscous energy dissipation is not 
taken into account as it will prove to be negligible with respect 
to the assumed flux from the crust or the energy release from 
nuclear burning. 

We consider the accretion of matter which is hydrogen- 
deficient, but has a large (99%) helium content. This allows us 
to study the stability of helium shell burning without the com- 
plication of the nuclear processes involving hydrogen. We will 
discuss the potential effect of the inclusion of hydrogen on our 
results. 

Relativistic effects are not considered during the simulation, 
since they are small in the thin helium burning shell. The effect 
of general relativity on observables, such as the luminosity and 
mass accretion rate, can be taken into account with appropriate 
gravitational redshift factors (e.g. lWoosley et al.ll2004) , reducing 
the values for the observables typically by several tens of per- 
cents. We do not apply these corrections to our results, as we are 
primarily interested in the relative effect of rotationally induced 
mixing on the stability of shell burning. 

As our models are one-dimensional, any asymmetries in 
the longitudinal and latitudinal directions are not resolved. For 



example, burning likely starts at one locatio n in the enve- 
lope, after which it spreads to other regions (B ildstenl 1 19951: 
Spitkovsky et al. 2002; Cooper & Narayan 2007). This idea is 
supp orted by the observation of so-called burst osc illations 
( Strohmaver etal.1 1 19961: IStrohmaver & Bildstenl 12006) . which 
are periodic variations in the light curve at the spin frequency 
of the neutron star that are observed during some type-I bursts. 
Also, Bildsten (1995) proposes a regime where the burning 
spreads relatively slowly across the surface on a similar time 
scale as the accretion of matter. We restrict ourselves to the case 
of stable helium burning, where flame spreading is not an issue. 

2.1. Grid 

The one-dimensional grid divides the star in shells with a cer- 
tain mass and is separated in two parts to facilitate matter ac- 
cretion while limiting numerical noise in the nuclear burning 
region. In the inner Lagrangian part each shell is assigned an 
absolute amount of mass, while in the outer so-called pseudo- 
Lagrangian part a shell is assigned an amount of mass relative to 
the total mass in the model. To minimize numerical diffusion, we 
choose these regions such that the nuclear burning takes place in 
the Lagrangian part. After each time step the grid is adjusted 
to ensure that any gradients in for instance temperature, energy 
generation rate or chemical composition are well resolved. 

Apart from the different grid points our model has the notion 
of a "core", which represents the inner part of the star, that is not 
modeled in detail and which sets the inner boundary conditions 
of our grid. The core has a mass, radius, rotation rate and a lumi- 
nosity. These quantities are kept at fixed values. During consec- 
utive time steps angular momentum is exchanged between the 
grid and the core, after which we return the angular momentum 
and velocity of the core to the initial values, effectively prevent- 
ing the core to spin up or down. We can safely make this approx- 
imation, since the amount of angular momentum that is accreted 
on the timescale of nuclear burning is negligible with respect to 
the angular momentum of the core. 

2.2. Accretion 

Mass accretion is implemented by increasing the total amount 
of mass and scaling the mass of the grid points in the pseudo- 
Lagrangian region. The mass fractions of the isotopes at each 
of these grid points are set to the weighted mean of the original 
composition of the grid point and the assumed composition of 
the accreted matter, where we weigh by, respectively, the original 
mass of the grid point and the accreted mass that is added to the 
grid point. 

The accreted matter originally carries angular momentum 
corresponding to orbital motion at the Keplerian angular veloc- 
ity at the surface of = VGMR~^, where M and R are the 
mass and radius of the star, respectively. We assume, however, 
that the angular momentum of the accreted material is quickly 
shared with a number of the outermost grid points. This region 
coincides with or is slightly smaller than the pseudo-Lagrangian 
region where the mass is added. The distribution of the angular 
momentum is done such that the specific angular momentum is 
constant in the outer region. After the accreted angular momen- 
tum is added to the model, the mixing of angular momentum 
resulting from the different viscous processes is applied to the 
model. As in the treatment of the accretion of mass, we make 
sure that the region where angular momentum is accreted lies 
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outside the region of interest, i.e. well above the helium burning 
shell. 



2.3. DifFusivity and viscosity 

Due to the accretion of angular momentum, the outer layers of 
the star spin faster than the core. Therefore a non-zero shear q is 
present 



dln^l 
dlnr ' 



(1) 



where r is the radial coordinate. The shear can give rise to hy- 
drodynamical instabilities. We consider the dynamical shear in- 
stability, the secular shear and the Solber^-H0i land in stability as 
well as Eddington-Sweet circulation (see, e.g. jHegeT et al. 2000 
for details of the implementation). Furthermore, we consider the 
presence of a rotationally induced magnetic field. 

A magnetic field is employed followin g the presc ription of 
the "Tayler- Spruit dynamo" presented by SpruitI (12002 V A small 
initial radial magnetic field Br winds up due to the shear to cre- 
ate a much stronger toroidal field . The field growth is limited 
by magnetohydrodynamic (MHD) instabilities. In this process 
Br is increased as well, such that once an equilibrium situation is 
reached, the initial field strength Br is unimportant. Equilibrium 
is reached on much shorter timescales than the mixing processes 
induced by the magnetic field. Therefore, in our model we em- 
ploy the equilibrium values for Br and B^ . 

Both the magnetic and non-magnetic hydrodynamic insta- 
bilities cause turbulence which allows for chemical mixing and 
the transport of angular momentum. Furthermore, angular mo- 
mentum is transported as well due to magnetic torques. We treat 
the chemical and angular momentum mixing as diffusion pro- 
cesses. Transfer of angular momentum is governed by the one- 
dimensional diffusion equation 



1 d 

pr^ dr 



da 



(2) 



with V the viscosity and p the density. Diffusion of the mass 
fraction 7 of a given isotope, e.g. ^He, is calculated via 



dY 
'dt 



1 d f 2dY\ 



(3) 



where D is the diffusivity. D and V are measures for the ef- 
ficiency of the respective diffusion processes. D is generally 
smaller than v, since it requires more work to exchange material 
than to exchange angular momentum. For the non-m agnetic in- 
stabilities we use V = 30D (IChabover&Zahnlll992h . This is not 
the case for the D and v that follow from the Tayler- Spruit mech- 
anism, where angular momentum is transported due to magnetic 
torques, such that D and v differ in their dependencies on q and 

The Tayler-Spruit formalism accounts for both composi- 
tional and a temperature gradients, which provide a stabilizing 
stratification against the magnetohydrodynamic instabilities that 
limit the growth of the magnetic field. The effective radial vis- 
cosity Vre is described as a function of Vgo and Vgi , the viscosities 
when thermal diffusion either can be ignored or dominates, re- 
spectively: 



VePVel 
VeO + Vei 



(4) 



where the factor < f{q) < 1 is introduced to ensure v^^van- 
ishes smoothly below a certain minimum shear (ISpruijl2002h . 



Veo = rW(f)' 



(5) 



with and A^r respectively the "compositional" and t hermal 
part of the buoyancy frequency and k the opacity (see I SpruitI 
I2OO2 for more details). Analogously, 'Spruit' ('2002') derive the 
following expression for the effective radial diffusivity Dr^: 



^eO^el 
^eO+^el 

- ^2 



m. 



(6) 



Analytic studies often use as approximation only Vgo, Dqq or Vgi , 
Del . Since we study thermonuclear burning in a thin shell, both 
thermal and compositional gradients are expected to be impor- 
tant in the shell. In our models we, therefore, use expressions (|4]) 
and (|6]) including both parts. 

While the Tayler-Spruit dynamo is often applied in the spher- 
ical approximation, w e use it here for a thin shell (see also 
iPiro & Bildstenll2QQ7l) . For the dynamo to work, the presence of 
a non-zero shear locally is sufficient. Since the resulting mag- 
netic field is predominantly toroidal, the angular momentum 
transport and the chemical mixing induced by the dynamo are 
limited to the radial region where the shear is present. 

In Section |4J1 we discuss some potential problems with the 
Tayler-Spruit mechanism. 

2.4. Nuclear network 

An extensive nuclear network consisting of 35 isotopes and 65 
nuclear reactions is employed to model the chemical evolution. 
The most important reaction sequence for our models is the triple 
alpha process, for which the net reaction is: 



34He 



.12 



C + y. 



(7) 



Furthermore, alpha capture reactions onto carbon and oxygen 
are considered: 



4He 
4He 



_16o. 



>i6o + 7 
Ne 



(8) 



We disable the energy output of the ^^C + ^^C and ^^C + ^^O 
reactions to prevent the occurrence of unstable carbon burning 
at larger depths in our models. Any heating that stable carbon 
may have provided can be modeled as an e xtra contribution to 
the heat flux from the crust (see Section [277]) . 

2.5. Viscous and compressional heating 

Viscous energy release due to the presence of shear heats the 
neutron star envelope. The energy generation rate per unit mass 
by viscous heating, eyisc. can be calculated as 



(9) 



We find that the relatively small shear present in our models, 
in the range of rotation rates and mass accretion rates we con- 
sider, leads to a negligible amount of viscous heating, as the to- 
tal energy generated per unit time is several orders of magnitude 
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smaller than the assumed flux from the crust or the energy gen - 
erated by stable helium burning (see also lPiro & Bildstenll2Q07l) . 
Therefore, we do not take viscous heating into account in our 
calculations. 

The continuous accretion of matter onto the neutron star 
compresses the matter in the envelope. This gives rise to com- 
pressional heating, which is calculated in our models. 

2.6. Initial model 

We start with a grid of approximately 600 grid points that rep- 
resent the outer layers of the neutron star from a column depth 
ofyc^ 1 gcm~^ down to j :^ 10^^ gcm~^, with a total mass of 
0^ 10^"^ g. As initial composition we take 50% ^^C and 50% ^^O. 
This acts as a buffer which is inert with respect to the triple-alpha 
reactions of helium burning. On this grid matter is accreted with 
a composition of 99% ^He and 1% other isotopes, based on the 
equilibrium mass fractions of the CNO cycle. 

We use a core with a mass of M = 1 .4M0 and a radius R = 
10 km for a canonical neutron star. The angular velocity of the 
core, is expressed in units of the Keplerian angular velocity 
at the surface: = 1.4 • 10^ radians s~^ We create models with 
different values for Q. 

The heat flux from the neutron star crust is modeled by set- 
ting the luminosity at the inner boundary of our grid to a fixed 
value Lcrust- 

2. 7. Flux from the crust 

As the ignition of the accreted helium depends strongly on the 
temperature and we do not consider heating by stable hydro- 
gen burning, it is influenced by Lcmst- The crust is a deeper 
layer of the neutron star starting at j lO^^gcm"^ and extend- 
ing down to the core where densities exceed the nuclear den- 
sity. When matter is accreted onto the neutron star, the crust 
is compressed. This gives rise to pycnonuclear burning and 
electron-capture reactions which heat both the neutron star core 
as well as the outer layers. Therefore, this heat flux is gener- 
ally assumed to scale with the accretion rate and a fixed energy 
generated per accreted nucleon is applied. P revious studies of- 
ten u s e a value close t o 0.1 M eV/nucleon (iHaensel & ZdunikI 
l2003h . [Gumming et al." (2006) use values in the range of 0.1- 
0.3 MeV/nucleon when studying the ignition conditions of long 
duration type I X-ray bu rsts. T he analysis of the superburst from 
4U 1608-52 (.Keek et al.l l2008') indicates that the temperature of 
the crust was likely much higher than predicted by current mod- 
els. Recently 'Gu pta et aD (2007) calculated the thermal structure 
of the crust, for the first time including electron captures into ex- 
cited states, finding a heat flux from the crust up to ten times 
highe r than obtained in previous studies (see also iGupta et aP 
l2QQ8h . It is clear that the value of the heat flux from the crust is 
still not well determined. Therefore, we choose to vary it inde- 
pendently from the mass accretion rate. 

Heat from helium burning in the envelope may be conducted 
into the core. From here it is either dissipated by neutrino emis- 
sion or conducted outwards, back to the envelope on a thermal 
time scale of the order of 10^ yr. This gives an extra contribu- 
tion to Lcrust, the size of which depends on the amount of heat 
from the envelope conducted inwards and on the efficiency of 
neutrino cooling of the core. The latter is determined fro m the 
equat ion of state of the core (see discussion in Gumming et al.l 
20061). which is not well constrained at this moment (see, e.g., 
Lattimer & Prakashll2007h . Therefore, a certain value for Lcmst 



can physically be obtained by different combinations of an out- 
ward flux from the core and the crust as well as an inward heat 
flux from the envelope. In our calculations we do not go into 
these details, but treat Lcmst as a free parameter. When we make a 
comparison to the physical situation, we will assume that crustal 
heating gives the domination contribution to Lcmst- 

By choosing a certain value for Lcmst, we fix the total heat 
flow into and out of the envelope to a net flux into the envelope. 
For studying unstable burning, a flash would give rise to a large 
heat flow out of the envelope and our boundary condition would 
not be valid for the duration of the flash. However, we consider 
only stable burning, which leads to a relatively lower heat flux 
that is constant in time and is compatible with our boundary con- 
dition. 

3. Results 

3.1. Turbulent mixing 

The accretion of angular momentum gives rise to shear, which 
causes turbulence and induces a magnetic field. These drive the 
transport of angular momentum and the mixing of the chemi- 
cal composition. In turn, the chemical and angular velocity pro- 
files influence the diffusivity and viscosity, until an equilibrium 
is reached. We present the properties of equilibrium models with 
different rotation rates of the core and different mass accretion 
rates. The heat flux from the crust is set to a value such that 
steady- state burning is achieved. 

Fig. [T] shows the rotational velocity profile for several values 
of the mass accretion rate. The viscous processes transport angu- 
lar momentum efficiently through the outer layers of the neutron 
star into the core, resulting in a small difference in between the 
outer and inner part of our model of ~ 10~^. The largest 

change in is at y lO^gcm"^. We will see that at this depth 
helium burning generates heat and changes the chemical com- 
position. In the figure we indicate the region of the model where 
the accreted angular momentum is added as described in the pre- 
vious section. Since the accreted angular momentum is quickly 
shared with this region, the angular momentum transport is not 
calculated in a physically correct way in that part of the model. 
Furthermore, at j < lO^gcm"^ lies the pseudo-Lagrangian re- 
gion where the accreted matter is added. Therefore, we will limit 
most plots to show the Lagrangian part, y > 10^ gcm~^. 

For higher accretion rates, angular momentum is added faster 
to the star, resulting in a higher shear (Figure[2]top), although the 
difference is less than an order of magnitude in the range of ac- 
cretion rates we consider. Since the accreted material originally 
rotates with the Keplerian velocity, models with a more slowly 
rotating core experience an up to 10^ times larger shear (Fig. [21 
bottom). The shear profiles show two broad "bumps": one in the 
region 10^ ^ J ^ lO^gcm"^ and another at j > lO^^gcm"^. We 
will see below that due to nuclear burning and diffusion there 
is a gradient in the chemical composition in these regions. This 
gradient has a stabilizing effect which reduces the angular mo- 
mentum transport, causing a somewhat larger shear. 

The presence of shear causes turbulence and induces a mag- 
netic field, that both drive diffusion and viscous processes. We 
obtain a magnetic field with a radial component 5^ :^ 10^ G and a 
toroidal component 5^ 10^^ G. We calculate D and V both due 
to the hydrodynamic instabilities and due to the (generation of) 
the magnetic field (Fig. [3l [4] and [5]). We find the latter to have the 
largest contribution to the diffusivity and viscosity. Of the non- 
magnetic instabilities, Eddington- Sweet circulation is the dom- 
inant process. At y lO^gcm"^ D and v due to this process 
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Figure 1. Rotational velocity ^1 as a function of column depth y 
for three values of the mass accretion rate and an angular veloc- 
ity of the core of 0. 1 ^1k- Accreted angular momentum is shared 
with the region to the left of the vertical line. 



peak, since at this depth the luminosity gradient is largest due to 
helium burning. Note that convective mixing plays no role in our 
models of steady state burning, because heat transport is radia- 
tive throughout our m odels, except at the onset of a flash, similar 
to what is found by Woos ley et al.l (l2004l) . 

We investigate how D and v depend on the rotation rate Q.. 
The contribution from the non-magnetic instabilities increases 
with increasing Q. (Fig.O. For the Eddin^ton-Swe et circulation 
the d iffusivit y and viscos ity are proportional to Q.^ (iKippenhahnI 
Il974j: [Heger et al. 2000). The diffusivity due to the magnetic 
field generation increases towards lower rotation rates as, ap- 
proximately, D oc Q,~^/^ (Fig. 131). The precise dependence is in- 
fluenced by buoyancy effects from temperature and composi- 
tional gradients (Sect. 12.31 Spruit 2002). The magnetic viscosity, 
which is induced by the magnetic torques instead of MHD in- 
stabilities, increases towards higher Q,, since it depends mostly 
on the rotation rate: D Qr^l^ (Fig. [5]). Interestingly, while 
the magnetic diffusivity dominates at most depths, for models 
with > 0.1 Eddington-Sweet circulation gives a contribu- 
tion to the diffusivity which is equal or larger at the depth of 
y 10^ gcm~^ where helium is burned. 

During our calculations we assume angular momentum is ac- 
creted with an efficiency of 100%. Due to viscous processes in 
the disk, the actual efficiency may be lower. We investigate the 
dependency of the shear, diffusivity and viscosity on the accre- 
tion efficiency by comparing a model with 100% to a model with 
0% efficiency. In the helium burning region, the difference is at 
most several tens of percents. For more realistic efficiencies the 
deviation from the 100% efficient model is likely even smaller. 
Therefore, we neglect this effect in our calculations. 

The viscosity allows for the transport of angular momentum, 
leading in the envelope to the angular velocity profile we pre- 
sented in Fig.[T] Similarly, diffusion mixes the chemical compo- 
sition of our model. Without diffusion, the chemical profile con- 
sists on the outside of the star (Fig. [6l middle, left-hand side) of 
accreted helium which burns aty 10^ gcm~^, forming a layer 
of mostly carbon. On the inside (right-hand side) the composi- 
tion is that of the initial model: 50% ^^C and 50% ^^O. When we 
include diffusion, ^He is mixed down to higher column depths , 
while some ^^C is mixed all the way to the surface (Fig. [6] bot- 
tom). Due to the increased helium abundance at j > 10^ gcm~^, 
the energy generation rate is higher there compared to the mod- 
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Figure 2. Comparison of shear q with column depth y. Top: 
models with different mass accretion rates and their rotation 
rate of the core fixed to ^1 = 0.1 ^^k- Note that the model with 
M = 0.6MEdd accreted for a longer time and, therefore, reaches 
larger depths. Bottom: models with different rotation rates of the 
core and a fixed accretion rate of M = O.SMEdd- The model with 
£1 = 0.9 suffers from numerical noise. 



els without diffusion (Fig. [6l top). For models with different ro- 
tation rates of the core, we see that those with higher diffusivity 
(see Fig. [4]) mix helium deeper inwards and consequently re- 
lease more energy due to thermonuclear burning at higher col- 
umn depths (Fig. [7]). 

In the steady state helium is burned typically at the depth 
where the accretion timescale equals the burning timescale. In 
the case of diffusive mixing of the chemical composition, fuel 
can be transported to deeper layers than the accretion column. 
In this case, helium is burned down to the depth where the burn- 
ing timescale equals the diffusion timescale. We determine as 
a function of column depth the burning and chemical diffusion 
timescales ^nuc and ^diff for ^He by calculating 



/AY 
^' AT 



(10) 



with Y the "^He mass fraction and AY the change in Y due to 
nuclear burning or chemical diffusion during one time step of 
our calculations (Fig. [8]). In the inner part of the burning zone, 
where helium is mixed in, the chemical diffusion and burning 
timescales differ by at most a few tens of percents, ensuring that 
steady- state burning can take place down to larger depths than 
without diffusion. 
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Figure 3. Non-magnetic diffusivity D and viscosity v due to 
Eddington- Sweet circulation as a function of column depth y for 
different values of the angular velocity of the core Q. and an ac- 
cretion rate of M = O.SMEdd- Note that v is larger than Z) by a 
fixed factor of 30. Especially the model with = 0.9 suffers 
from numerical noise. 
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Figure 5. Magnetic viscosity v as a function of column depth y 
for different values of the angular velocity of the core Q,. The 
model with = 0.9 suffers from numerical noise. 
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Figure 4. Magnetic diffusivity Z) as a function of column depth 
y for different values of the angular velocity of the core and an 
accretion rate of M = O.SMEdd- A minimum in D is apparent at 
^1 0. 1 ^Ik- The model with Q, = 0.9 Q.^ suffers from numerical 
noise. 



3.2. Effect of mixing on stability 

The temperature in the burning shell has an important influence 
on the stability of the nuclear burning. Typically, burning is sta- 
ble for higher temperatures of the neutron star envelope, which 
depends highly on the heat flux from the crust. We express this 
heat flux as a luminosity Lcmst and investigate for which values 
burning is stable by starting at a high value of Lcmst and subse- 
quently lowering it at a rate which is slow with respect to the 
thermal timescale of the inner part of our model. For a certain 
value of Lcmst nuclear burning no longer proceeds stably and a 
flash occurs. We refer to this as the critical luminosity Lent- At 
that point we stop the calculation. 

We perform this procedure for two models with £1 = 0.1 
and M = O.SMgdd- Both models take viscous processes into 
account, such that the angular velocity profiles are similar. 
However, we enable chemical diffusion only for one model, 
while the other is calculated without chemical diffusion. For the 
model excluding chemical mixing, we start with stable burn- 
ing (Fig. [9] top). Lowering Lcmst, helium shell burning exhibits 
oscillations with increasing amplitude. These oscillations are 



bJO 



1015 

1014 
1013 
IOI2 

lOii 
lOlo 

109 

108 



no diffusion 

~ ~ diffusion 



1 \ \ \ r 



J \ \ \ \ L 





10^ 
10-1 

10-2 

10-3 
10- 



■ 


1 1 1 1 1 1 1 1 
diffusion / \ 1 


\ 1 1 ■ 



10° loi 10^ io3 10^ 10^ 10^ 10'^ 10^ 10^ 101010111012 
y (g 



Figure 6. Nuclear energy generation rate Cnuc and isotope mass 
fractions as a function of column depth y for two models with 
M = O.SMEdd and Q. = 0.1 ^^k- One model is calculated includ- 
ing chemical diffusion, while the other is not. 



well resolved in time by our calculations. For Lcmst = 2.448 • 
lO^^ergs"^ the oscillation amplitude increases rapidly, until a 
flash occurs. We regard this as the onset of unstable burning. 
It is possible that what we regard as the onset of instability is 
actually just an oscillation with a large amplitude, which could 
be followed by a series of oscillations before the unstable situa- 
tion is reached. Nevertheless, observing how fast the amplitude 
of the oscillations increases in a small interval of Lcmst, we ex- 
pect that our measure of Lcrit is accurate within 0.2%. The period 
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Figure?. Nuclear energy generation rate £nuc as a function of 
column depth y for different values of the angular velocity of the 
core Q.. Models with higher diffusivity (see Fig. [4]) have more He 
burning at higher column depths. The model with = 0.9 
suffers from numerical noise. 
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Figure 8. Timescales for helium burning, chemical diffusion and 
accretion as a function of column depth y for Q. = 0.1 and 
M = 0.3MEdd. 



of the pulses increases, from approximately 6 • lO^s at Lcmst = 
2.448 • 10^^ ergs ^ to 9 • lO^s at Lcrust = 2.444 • lO^^ergs \just 
before the final two oscillations. 

For the model including diffusion, the behavior is similar 
(Fig. [9] bottom). The amplitude of the oscillations, however, is 
large for a broader range of Lcmst and the behavior is less reg- 
ular. During intervals of 10^ — 10^ s the amplitude slowly in- 
creases and subsequently decreases either slowly or fast. The 
cause of the irregularities is likely the increased complexity of 
our model when turbulent mixing is included, giving rise to more 
"noise". Despite the irregular behavior, we can still discern a 
general trend towards oscillations with larger amplitudes while 
we lower Lcmst- The period of the oscillations is ~ 1.2- 10^ sand 
does not change significantly. The value of Lcrit we obtain for 
this model is over a factor two smaller than for the model without 
diffusion: Lcrit = 1.180- lO^^ergs"^. This shows the substantial 
stabilizing effect of turbulent mixing. The uncertainty in Lcrit is 
likely larger than for the model without diffusion due to the in- 
creased noise. Note that for the chosen angular velocity of the 
core, = 0.1 ^1k, the diffusivity is at its minimum value. For 
another choice of ^1, the effect is likely to be even larger. 

We investigate the dependence of Lcrit on by calculating 
two series of models, including and excluding chemical dif- 
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Figure 9. Nuclear luminosity Lnuc due to helium burning as a 
function of time, while lowering Lcrust at a constant rate. Shown 
are the last 10^ s before the onset of unstable burning. The in- 
sets zoom in on the last 10^ s. Top: model without diffusion for 
£l = 0.1 ^^K, M = O.SMEdd- While burning is initially stable, we 
see an increase in the amplitude of the oscillations, until finally 
burning becomes unstable. Bottom: similar model including dif- 
fusion. The amplitude of oscillatory burning is much larger and 
less regular than for the model without diffusion. 



fusion. Due to the increased noise, the models which include 
diffusion are computationally more expensive. Therefore, we 
limit ourselves to two models with relatively low rotation rates, 
= O.OS^^K, O.l^^K, and two models with high rotation rates: 
= 0.8^1k, 0.9^1k- We find that Lcrit is lower by more than a 
factor of 2 for all values of we consider (Fig.[TO|. Furthermore, 
even for the models without diffusion we see a decrease in Lcrit 
of up to 12% when going to higher rotation rates. This stabilizing 
effect is due to the centrifugal force F^f ^ Q^^, which lowers the 
effective gravity, increasing the efficiency of stabilizing effects. 
Note that due to the one-dimensional approximation we use, this 
is likely an underestimation (see Sect. 14.21) . 

We noted earlier that diffusivity due to Eddington- Sweet cir- 
culation can be comparable to or even exceed the diffusivity 
due to magnetic fields. To investigate the effect of non-magnetic 
instabilities alone, we created a model with = O.OS^^k and 
M = O.SMEdd including only angular momentum transport and 
chemical diffusion due to non-magnetic instabilities. Angular 
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Figure 10. The heat flux from the crust at the transition of stable 
to unstable burning Lcrit as a function of the angular velocity 
of the core in units of the Keplerian velocity at the surface for 
models with M = O.SMEdd- Each dot represents one model. Two 
series of models are calculated: without (upper) and with chemi- 
cal diffusion (lower). One model is calculated without the effect 
of magnetic fields (middle). Here the secular shear instability 
drives the chemical mixing. On the right-hand vertical axis is 
indicated the heat generated in the crust per accreted nucleon, 
Gcrust, for the given accretion rate. The dotted line indicates 
Gcrust = 0.1 Me V nucleon" \ which is used in many studies. 



momentum transport is much less efficient, leading to a large 
shear. The secular shear instability drives the diffusivity and vis- 
cosity. For this model we find Lcrit = 1.948 • 10^^ erg s-^ which 
is 23% lower than what we find without diffusion (Fig.fTOb. 

When we change Lcmst it takes a thermal timescale for our 
model to adjust. Since we lower Lcmst continuously, this intro- 
duces a systematic error in the value of Lcrit we obtain, as the 
model is not fully adjusted to that value of Lcrust when we stop 
the calculation. To estimate the systematical error we take one 
model and lower Lcrust at different rates, all slow with respect to 
the thermal timescale of the model, and find that the difference 
in Lcrit and, ergo, the systematic error is about 4%. 

3.3. Mass accretion rate dependence 

The dependence of Lcrit on the mass accretion rate is investigated 
by calculating two series of models with and without chemical 
diffusion, in a similar way as we studied the dependence on Q,. 
The angular velocity of the core is set to O.IQ,^. We find that 
Lcrit increases towards higher mass accretion rates (Fig.[TT]). For 
the models without chemical diffusion, the accretion rate above 
which the luminosity of the crust exceeds Lcrit as expected from 
a crustal energy release of Qcmst = 0.1 Me V nucleon" \ which is 
used in many studies (see Sect. l2.7b , l ies at MEdd. wh ich is lower 
than found in previous studies (e.g., Bildsten"l998'). Including 
chemical diffusion, the transition takes place at lower rate: ap- 
proximately 0.4MEdd- 

To understand the behavior of increasing Lcrit for larger val- 
ues of M, we follow the stability analysis of therm onuclear 
burning in a thin s hell by Giannone & Weigert (*1967) (see also 
lYoon etal.|[2004bh . These authors derive an expression for the 
inverse of the growth time scale of temperature perturbations, 
'^pert' considering both thermonuclear burning and radiative 
cooling (equation 16 in Giannone & Weigert 1967). In this ex- 
pression the luminosity from below the thin shell is ignored. 
We can include Lcrust by adding it to the luminosity of the 
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Figure 11. The heat flux from the crust at the transition of sta- 
ble to unstable burning Lcrit as a function of the accretion rate 
M for models with = 0.1 ^Ik- Each dot represents one model. 
Two series of models are calculated: without (upper) and with 
chemical diffusion (lower). The dotted line indicates the lumi- 
nosity from the crust as a function of M for an energy release of 
Gcrust = 0.1 Me V nucleon" \ which is used in many studies. 



shell, which results in an extra t erm in the numerator of th e 
right-hand side of equation (16) in lGiannone & Weiger^ (1 19671) : 
TpgJ^ oc —Lcrust /AM. Here, AM is the mass of the shell, which for 
steady-state burning is proportional to the mass accretion rate. 
Increasing M leads to a lower perturbation growth time scale. 
Therefore, to keep the same perturbation growth time scale, a 
larger value of Lcrust is required at larger M. This simple picture 
suggests that Lcrit depends linearly on M, while Fig. [TT] suggests 
for the series of models without chemical diffusion that the re- 
lation deviates from a linear dependence. In fact, for larger M 
also the temperature and density in the shell are larger, which 
influences the perturbation growth time scale as well. 



4. Discussion 

We performed detailed calculations of thin shell burning of ac- 
creted helium on a neutron star. Similar to what was found by 
lYoon et aD (l2004al) for white dwarfs, also for neutron stars rota- 
tio nal mixing l eads to more stable burning. As was demonstrated 
by iPiro & B ildsten (2007), rotationally induced magnetic fields 
provide in mo st cases the domina nt contribution to the diffusivity 
and viscosity. Yoon et all (l2004al) considered only non-magnetic 
hydrodynamic instabilities. In those models the transport of an- 
gular momentum is less efficient, resulting in a much larger 
shear. The secular shear instability is in that case the driver of 
turbulent mixing. We take magnetic fiel ds into account through 
the Tayler- Spruit dynamo tepruit|[200 2), which causes efficient 
angular momentum transport and turbulent mixing. Eddington- 
Sweet circulation also plays an important role, especially in the 
burning zone, at the peak of the energy generation rate, where 
there is a steep gradient in the luminosity. 

In our models with diffusion the geometrical thickness of 
the burning shell is larger (Fig. [6|. For instance, for a model 
with = 0.1 and M = O.STl^Edd the thickness of the burning 
shell is 6.7m excluding and 8.7m including chemical diffusion, 
where we define the shell such that at the boundary the energy 
generation rate is 10~^ times the peak value. This results in an 
increased stability of thermonuclear burning. 

Due to the increased stability, the critical value of the heat 
flux for the crust Lcrit at the transition of stable and unstable burn- 
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ing is lower for models that include chemical diffusion. Even for 
the rotation rate which gives the minimum diffusivity, nuclear 
burning proceeds stably down to values of Lcrit that are a factor of 
two lower than for the non-diffusion case. In these calculations 
we vary the flux from the crust independently from the mass ac- 
cretion rate. For several values of M and we determine Lcrit- 
Observationally, this boundary is placed at M :^ 0. IMEdd- At this 
accretion rate, we find Lcrit = 0.508 • lO^^ergs"^ for a model 
with = 0.1 ^1k- Assuming the heat flux from the crust is the 
main source of this luminosity, this corresponds to an emergent 
heat from the crust of 0.27MeVnucleon~^ This is almost three 
times the value of O.lMeVnucleon"^ used in many previous 
studies. Nevertheless, it is in agreement wi th the increa s ed he at 
deposition in the crust from calculations by lGupta et al.l (l2007h . 

It is especially interesting to note that we obtain this result 
already at the v a lue fo r for which the diffusivity is minimal. 
iPiro & BildstenI (l2007l) find this behavior only at low rotation 
rates: < 0.1 ^^k- The Taylor-Spruit mechanism considers sta- 
bilizing effects due to temperature and compositional stratifica- 
tion. Piro & BildstenI ([2007) simplify their models by taking into 
account only the former when determining the recurrence time 
of bursts and the condition for stable burning. This allows for a 
simple dependence of D on of D oc However, by in- 

cluding both stabilizing effects the dependence of D on ^1 and q 
is more complex. In our calculations we find that the spread in 
the diffusivity in the considered range of rotation rates is at most 
one order of magnitude. Furthermore, we find that at high angu- 
lar velocities Eddington- Sweet circulations drive mixing just as 
efficiently as the magnetic instabilities at low rotation rates. This 
means that for all rotation rates the stability of shell burning is 
increased by chemical mixing due to rotation and rotationally 
induced magnetic fields. 

4.1. Mixing efficiencies 

In our models we find that the magnetic instabilities and torques 
provide the most efficient mixing of the chemical composition 
and angular momentum. The only exception is found at higher 
rotation rates, where the Eddington- Sweet circulations provide 
an important contribution to the diffusivity in the helium burning 
shell. Note that we employ a one-dimensional approxima tion of 
the di ffusivity due to Eddington- Sweet circulations ( Kippe nhahnl 
11974 . At higher rotation rates two-dimensional effects may be 
important that cannot be tested with our model, but that may 
influence the efficiency of mixing due to these circulations. 

Several studies need t he Tayler-Spruit d ynamo to repro- 
duce observed quantities: Heger et alT (l2005h model the pre- 
supernova evolution of massive stars and use the angular mo- 
mentum transport due to rotationally i nduced magnetic fields 
to explain observed pulsar frequencies. ISuijs et aLr(l2008 l) per- 
formed similar calculations for white dwarf progenitors and find 
that angular momentum transport by magnetic torques is needed 
to explain the observed range of white dwarf rotation rates. 

Ho wever, the ni e chanis m of the Tayler-Spruit dynamo is un- 
certain. IZahn et aP (l2007l ) performed calculations with three- 
dimensional models which do not reproduce the dynamo effect. 
Therefore, it is possible that this mechanism ha s a lower ef- 
ficien cy or is absent (see, e.g., the discussion in iMaeder et al.l 
I2008h .' In that case, non-magnetic instabilities drive chemical 
mixing and angular momentum transport. These are less effi- 
cient by many orders of magnitude. This leads to a large shear, 
causing the secular shear instability to provid e the domina nt con- 
tribution to D and v, similarly to what .Yoon et al.l (l2004a|) found 



in white dwarf models which did not include a rotationally in- 
duced magnetic field. We performed calculations on one model 
without a rotationally induced magnetic field and we find a value 
of Lcrit that is 23% lower than for similar models that do not take 
chemical diffusion into account (Fig.[TOl). This model has a ro- 
tation rate of O.OS^^k- Models with a slower rotating core will 
have a higher shear, which may result in a lower value of Lcrit • 
Therefore, even without magnetic fields, the effect of rotational 
mixing on the stability of helium shell burning is considerable. 

4.2. Centrifugal force 

Another stabilizin g effec t of rotation is provided by the cen- 
trifugal force. Yoo n et all (|2004a) found that for shell burning in 
white dwarfs, the centrifugal force lowers the density and degen- 
eracy in the shell, which leads to increased stability. We created 
a series of models without chemical mixing and found that for 
increasing Q,, Lent is lower by up to 12%. However, this likely is 
an underestimation. Our implementation of the centrifugal force 
is a one-dimensional approximation. The centrifugal force de- 
pends on latitude, being largest at the equator and vanishing at 
the poles. Our calculations use a value which is averaged over 
latitude (e.g. Heger et al. 2000). In the case of rigidly rotating 
polytropes this approxi mation has been f ound accurate for rota- 
tion rates up to 0.6^^k ( Yoon et al.l2004al and references therein). 
Therefore, we limit the effect of the centrifugal force above 60% 
of the critical velocity. In actuality, the centrifugal force may re- 
duce Lcrit even further for fast spinning neutron stars. 

4.3. Marginally stable burning 

iRevnivtsev etaP (l200lh observe from the LMXBs 4U 1636-53, 
4U 1608-52 and Aql X-1 millihertz quasi-periodic oscillations 
(mHz QPOs) and speculate that this phenomenon may be related 
to a special mode of nuclear burning. This behavior is thought to 
occur at the boundary of stable and unstable burning, since it is 
observed in a narrow range of accretion rates, where at lower ac- 
cretion rates type I X-ray bursts were observed while they were 
absent at higher accretion rates. The oscillatory burning mode at 
this bo undary was expected from one- zone models by Pac zynskj 
(Il983h as "marginally stable burning". iHeger et al.l (|2007|) repro- 
duce this behavior in multi-zone numerical models of the neu- 
tron star envelope and derive the period of the oscillations Pqsc to 
be the geometric mean of the accretion and thermal timescales: 

^os c A/^acc^therm- 

lAltamirano et al.l (|2008|) observe frequency drifts in mHz 
QPOs from 4U 1636-53. The mHz QPOs are observable for ap- 
proximately lO^s. During this time the frequency drops by sev- 
eral tens of percents. At a certain point a type I X-ray burst is 
observed, after which no oscillations are detected for some time. 
The change in period is not found to be directly related to vari- 
ations of the X-ray flux on short timescales. Therefore the fre- 
quency of the oscillations does not appear to be sensitive t o vari- 
ations in the accretion rate. Furthermore, iMolkov et al.l (l2005l) 
observed a long type I X-ray burst from SLX 1735-269 which 
exhibits ~ 10 mHz oscillations during the burst decay. Also in 
this case a drift to lower frequencies takes place during the part 
of the X-ray light curve where the black body temperature is 
observed to decrease. 

Also in our models we observe oscillations at the boundary 
of stable and unstable burning. This is most pronounced in the 
models without chemical diffusion. When we lower Lcrust, effec- 
tively cooling our burning layer, the frequency of the oscillations 
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decreases by several tens of percents until a flash occurs. The fre- 
quency drift takes place over a small range of Lcmst, just before 
the onset of a flash. The amplitude of the oscillations is larger in 
this range, making the osc illations more eas ily detectable. This 
supports the suggestion by lAltamirano et al.l (.2008) that the fre- 
quency drifts are due to cooling of deeper layers and could ex- 
plain why it is not sensitive to short-term variations of the accre- 
tion rate. The cooling could, for instance, be caused by the slow 
release of energy from a deeper layer that was heated up during 
an X-ray burst. 

Note that the models that include diffusion display oscilla- 
tions as well, but the behavior is much more irregular, as the 
amplitude of the oscillations changes in a seemingly erratic way. 
This makes it difficult to draw conclusions from the oscillatory 
behavior. We suspect that numerical noise is the source of this 
behavior. A different implementation of the mixing processes in 
the model may be required to improve this. 



4.4. Further consequences of mixing 

Apart from the increased stability of helium burning, chemical 
mixing can give rise to other effects. Superbursts, which are en- 
ergetic type I X-ray bursts that last up to a day, are thought 
to be the result of a carbon flash ( Strohmayer & Brown 2002; 
Gumming & Bildsteni l2001l) . Fifteen superbursts have been ob- 
served from ten sources, most of which have accretion rate s M > 
Q.lMFdd fc ulkers 2Q0l lin 't Zand et al.ll2b 04: Remillard^etalj 
l2005HKuulkersii2005>) . As noted by'Piro & Bildsten (2007), tur- 



bulent mixing could allow for stable burning of helium at these 
accretion rates, building up a thick carbon layer. However, all 
superbursters exhibit normal type I bursts as well. In the pres- 
ence of hydrogen, the rp-process creates heavy isotopes, while 
lowering the carbon abundance. It has proven difficult to explain 
how carbon can surviv e many normal bursts to fuel a superburst 
(^e.g. JKeek et"ani2008h . 

While helium is mixed down, we find that carbon is mixed up 
to smaller depths (Fig. [6]). For our helium- accreting model this 
is of little consequence, but many neutron stars in LMXBs ac- 
crete hydrogen-rich material. Hydrogen burns via the hot-CNO 
cycle. The energy generation rate of this process depends on the 
mass fractions of the CNO seed nuclei, which are increased if 
carbon is mixed into the hydrogen/helium layer. This may lower 
the hydrogen mass fraction in the burning shell. If subsequently 
a hydrogen/helium flash occurs, there is less hydrogen present 
for the rp-process to destroy the carbon created by helium fu- 
sion. This could allow for the creation of a thick carbon layer to 
fuel a superburst. Further calculations are needed to assess the 
magnitude of this effect. 

Our models show that carbon is mixed all the way to the 
neutron star photosphere aty = 1 gcm~^. In case hydrogen-rich 
material is accreted, the ashes of hydrogen and helium burning 
can contain heavy isotopes, e.g. iron. If mixing increases the 
mass fractions of these isotopes near the surface, this would al- 
low for the detection of absorption lines from the neutron star. 
Measuring the gravitational redshift of the absorption line pro- 
vides constraints on the neutron star equation of state, such 
as has been shown with the tentative result of ICottam et al.l 
("2002) (D^ 200^). This is analogous to the calculations of 
.Weinberg et al. ( 2006), who predict that convection during par- 
ticularly luminous X-ray bursts can mix heavy elements to the 
photosphere. 



5. Conclusion 

A long-standing problem for the theory of thermonuclear burn- 
ing in the envelope of accreting neutron stars is the accretion 
rate at which the burning changes from unstable to stable. While 
observations point to M O.lMEdd, the models find M M^dd 
when hydrogen-rich material is accreted and M lOMEdd when 
i t is hydrogen-deficien t. Recent analytical and numerical models 
(Piro ABildien''2007") suggest that turbulent mixing due to ro- 
tationally induced magnetic fields leads to an increased stability 
of the thin shell burning on neutron stars, similar to the effect 
of mixin g due to hydrodynamic instabilities in the white dwarf 
case ( Yoo n et al.| [2004a). 

We perform for the first time detailed calculations of nuclear 
burning of accreted helium in the neutron star envelope using 
an extensive numerical model including both magnetic and non- 
magnetic hydrodynamic instabilities. In most cases the mag- 
netic instabilities and torques provide the largest contribution to 
the diffusivity and the viscosity, respectively, while Eddington- 
Sweet circulations become important at high rotation rates. We 
create models with an accretion rate of M = O.SMEdd and dif- 
ferent rotation rates of the core. We do not couple the heat flux 
from the crust to the mass accretion rate, but leave it as a free pa- 
rameter, finding the threshold value between stable and unstable 
burning. This value is lower by at least a factor of two for models 
including diffusion, in comparison with models without chemi- 
cal mixing. If we ignore the rotationally induced magnetic fields, 
mixing due to the secular shear instability can lower the critical 
heat flux by several tens of percents. For larger rotation rates the 
centrifugal force allows for stable burning at crustal heat fluxes 
that are up to 12% smaller in our models. Note, however, that this 
is likely an underestimation, because of the employed approxi- 
mation of the centrifugal force in our one-dimensional model. 

These effects show that the stability of helium shell burn- 
ing is increased due to rotation, such that stable burning can 
take place at sub-Eddington accretion rates. Magnetic field gen- 
eration through the Tayler- Spruit dynamo yields the strongest 
effect. However, even without magnetic fields, hydrodynamic 
instabilities and the centrifugal force increase stability substan- 
tially. We calculate models including chemical diffusion for sev- 
eral values of the accretion rate and find that at approximately 
40% of the Eddington-limited rate the critical heat flux from 
the crust equals the value inferred from an energy release of 
Gcrust = 0.1MeVnucleon~\ which is used in many previous 
studies. Therefore, to explain the observed transition at O.lMEdd 
an increased energ y release is required, as found recently by 
iGupta et aP (l2007h . Combined with this larger energy release, 
rotational mixing and rotationally induced magnetic fields may 
explain the transition from unstable to stable burning at the ob- 
served accretion rate. A better understanding of the mixing effi- 
ciencies is required to further constrain the importance of these 
mixing processes. 

Our models exhibit oscillations close to the transition of sta- 
ble to unstab le burning, \yhich h ave been related to the observed 
mHz QPOs (iHeger et al.ll2007h . When we lower the heat flux 
from the crust at a constant rate, we find that the frequency of 
the oscillations is reduced by u p to several tens of percents. This 
supports the suggestion by Altamirano et al.l (12008) that the fre- 
quency drifts observed in mHz QPOs are due to cooling. 
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